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^""^ Abstract 

o 

CN| We study a reaction-diffusion system of partial differential equations, wlrich can be taken 

O to be a basic model for criminal activity, first introduced in 3J . We show that the assumption 

1^ of a populations natural tendency towards crime significantly changes the long-time behavior 

tin of criminal activity patterns. Under the right assumptions on these natural tendencies we 

QQ first show that there exists traveling wave solutions connecting zones with no criminal activity 

T—i and zones with high criminal activity, known as hotspots. This corresponds to an invasion of 

criminal activity onto all space. Second, we study the problem of preventing such invasions 

p ^ by employing a finite number of resources that reduce the payoff committing a crime in a 
finite region. We make the concept of wave propagation mathematically rigorous in this 

. situation by proving the existence of entire solutions that approach traveling waves as time 

i^H approaches negative infinity. Furthermore, we characterize the minimum amount of resources 

Cd necessary to prevent the invasion in the case when prevention is possible. Finally, we apply 

^ our theory to what is commonly known as the gap problem in the excitable media literature, 

I— I proving existing conjectures in the literature. 

^ 1 Introduction 
cn 

The work of Short and collaborators on the mathematical modeling of criminal behavior in |19j 
has generated much interest, not only in the mathematical community, but in other academic 
C<l and non-academic spheres. In [TU] the authors introduced a parabolic system of PDEs to study 

the dynamics of crime hotspots, which are spatial and temporal areas of high density of crime. 
I This work has been followed by a vast amount of research geared towards using mathematical 

j> tools in an effort to gain some understanding of the phenomena of crime, see for example [HI [15]. 

This includes the work of Berestycki and Nadal in [3], where they introduce a more general class 
of reaction-diffusion system of equations. In this work the authors introduce the concept of a 
warm-spot, which are areas with some criminal activity, but not enough to be classified as a 
hotspot. The use of PDE systems to model this phenomenon provides a different, interesting, 
and insightful perspective. Furthermore, a PDE model allows for additional mathematical tools 
not available in more commonly used discrete models. The focus of this work is the study a 
class of reaction diffusion systems of partial differential equations introduced in [3], as means 
of continuing to explore criminal behavior. In particular, we are interested in the role that 
population's natural tendencies towards criminal activity affect the long time patterns. Our 



X 



'corresponding author: nrodriguez@math.stanford.edu 



1 



starting point is the system introduced in [3]: 

st{x, t) = As{x, t) — s{x, t) + Sf){x) + {p{x) — c{x, t))u{x, t) (la) 

ut{x,t) = —{A{s)-u{x,t)) (lb) 

, , 1 / , , / u(x,t)p(x) \\ 

c,(x, t) = - ^-c(x, t) + CO + r?,co [j^^^^^;^^^ - 1 J J > (Ic) 

defined for x G M" and t > 0. The unknowns are: the moving average of crime, u{x,t), 
the population's propensity to commit a crime, s{x,t), and the cost of committing a crime, 
c{x, t). This model, as the model studied in [19], is based on two fundamental assumptions from 
criminology theory: routine activity theory p] and the repeat and near-repeat victimization effect. 
Routine activity theory states that the two essential elements for a crime to occur are a motivated 
agent and an opportunity, and neglects the other factors. This has the additional benefit of 
making the problem mathematically tractable. The repeat and near-repeat victimization effect, 
which has been observed in real crime data, is the effect that crime in an area leads to more 
crime [T21I2T] . 

We now briefly explain how ([T]) incorporates the above assumptions, see [3J for more details. 



Let us start with (lb). The most fundamental assumption is that the moving average of crime 
u{x,t) and the propensity to commit crime s{x,t) positively influence each other. If s{x,t) is 
non-positive then it does not affect the crime average but if s{x, t) becomes positive, its influence 
on the criminal activity increases. The fraction of the population at a given place and time that 
will commit a crime is given by 

= I 1 - e-^^ s > 0. (2) 
Here /3 > is a given constant. When s{x,t) = the crime will disappear on the time scale Tu, 



and the balance of these two effects gives ( lb ) 



The propensity to commit a crime, s{x,t), evolves proportionally to the amount of crime: it 
is observed in the crime data that crime leads to more crime, this is the 'repeat and near-repeat 
victimization effect'. Hence, we assume that s{x,t) increases proportionally to u{x,t) with a 
rate that is the difference of two effects: p{x) measures the payoff of committing a crime and 
c(x, t) is its cost. Note that p{x) is naturally spatially heterogeneous as different neighborhoods 
have different payoff for a successful crime. A good example to keep in mind is residential 
burglaries, where the net payoff of a successful burglary is higher in wealthier neighborhoods. 



In addition, the diffusion term in (la) takes into account the non-local effect that can also 



be seen on the example of residential burglaries. The constant si^ in (la) is the equilibrium 
value, or a measure of the population innate tendency towards criminal activity. If s;, < 
the population is assumed to have a natural anti- crime tendency, Sh = ^ assumes a natural 
indifference towards criminal activity, and s;, > assumes an natural tendency towards criminal 
activity. The function Sh{x) corresponds to the base attractiveness value Ao{x) in the model 
introduced in [19J. To the authors' knowledge there have been no studies of the important effect 
that the base attractiveness value has on the long time behavior of the solutions. Here we show 
that Sf, plays a significant role in the long term behavior of criminal activity patterns. 



Equation ( Ic ) models what is referred to in [3] as an adaptive cost, assuming that resources 
are allocated at high crime rate areas. This is known as hotspot policing. As a first step, we 
study a model that does not consider the dynamics of the cost ([s]) 

st{x, t) = As{x, t) — s(x, t) + Sfo + a{x)u{x, t) (3a) 
ut{x,t) = A{s) — u{x,t). (3b) 

Furthermore, we assume that Sb is spatially homogeneous (sf,(x) = Sb). We plan to address the 
full model in the near future. 
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We consider initial conditions that only depend on one spatial variable, xi, 



s(x,0) = so(xi) (4a) 
u(x,0) = uo(xi). (4b) 

When a < the system decays to equilibrium, and the dynamics is not particularly interesting, 
so we assume that a{x) > 0, which is also the more realistic case sociologically. 
A steady state solution satisfies 

u = A(s) and Sxx + G{s) = (5) 

where 

G{s) = -s + Sb + aA{s). (6) 

We see that indeed the sign of Sb determines the number of solutions to ^ and thus affects 
the long time behavior of solutions. Equation ^ provides a natural connection between the 
parabolic system ([s]) and the single parabolic equation 

St = Sxx + fis). (7) 

The parabolic equation ([T]) has been the subject of much study since its introduction by Fisher 
in [9j as a model for the spread of advantageous genetic traits in a population. Since then ^ 
has been used to model numerous other phenomena, from flame propagation [I] to nerve pulse 
propagation [16]. 

We consider three aspects of ([3]): first, the invasion of high criminal activity into areas with 
low initial criminal activity ~ this is related to the existence of traveling wave solutions. Second, 
we study the prevention of such invasions through what is commonly known, in the excitable 
media literature, as the gap problem [13]. In this formulation, limited resources are employed in 
a finite region in order to reduce the strength of reaction term (which fosters propagation). We 
prove that, in the right parameter regime, it is possible to prevent the propagation of criminal 
activity. Furthermore, we characterize the minimal length of the interval, where the reaction 
term is replaced by decay, required to prevent the crime wave propagation. We also prove the 
existence of a unique entire in time solution that approaches a traveling wave solution in the 
limit t —7- —CO. This result rigorously defines the notion of wave propagation in the heterogenous 
problem. Finally, we address the issue of splitting resources and prove that splitting resources 
is never beneficial. 

Previous work on preventing propagation of traveling waves in excitable media includes, for 
instance, applications in physiology see |17| and in chemical reaction theory [S]. Of particular 
relevance is the work of Lewis and Keener [13], who consider the gap problem for the single 
parabolic bistable equation ([T]) in M. Through a bifurcation analysis and a geometric approach 
they prove the existence of a unique critical length that stops the propagation of the wave. We 
generalize the results of [13] to a system of equations using different methods that seem to us 
to be more adaptable to other problems. In addition, we prove that there is a unique entire in 
time solution of the gap problem that converges to a traveling wave as t — )■ — oo. This result 
also applies to ([T]). 

Outline: We describe the main results in Section [2] The bistable system is treated in 
Section [3] and the monostable system is treated in Section |4j We devote Section [5] to the proof 
of existence of an entire solution to ([7]), which approaches a traveling wave solution as time 
approaches negative infinity. We conclude with a discussion in Section [6] 
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ematical Sciences DMS-1103765. Also, part of this work was carried out while Henri Berestycki 
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2 Main results 

2.1 Definitions, Notation, and Preliminaries 

In order to deal with non-negativity of solutions, we rewrite ([s]) in terms of s = s + si,- 

st{x, t) = As{x, t) — s{x, t) + aL{x)u{x, t) (8a) 
ut{x, t) = g{s) - u{x, t), (8b) 

where, Sh G M, g(s) = A(s — Sfe) and 

, . (a X < and x > L 
= I < X < L, 

for L > and a > 0. The case L = corresponds to the problem on a homogeneous medium 
and the case when L > corresponds to that in a heterogeneous medium. The steady state 
version of ([S]) is 

u = g{s) (10a) 
Sxx = -/(s), (10b) 

with 

f{s) = -s + ag{s). (11) 

Here we have replaced s with s for notational simplicity. Note that L > corresponds to the 
gap problem Q with /(s) replaced by —s in an interval of length L. 

The function /(s) has either one, two, or three zeros depending on the parameters Sh,a,/3. 
In the bistable case f has three zeros, and we assume without loss of generality that Sb, /3, and 
Q are such that there exists a unique a G (0, 1) so that 

/(O) = /(a) = /(I) = 0, < S6 < a < 1, and /'(O) < 0, f{a) > 0, /(I) < 0. (12) 

In this case, the steady state (0,0) represents zero criminal activity, {a,g{a)) is a warm-spot, 
and {l,g{l)) is a hotspot. The states (0,0) and {l,g{l)) are stable, and {a,g{a)) is unstable. 
Thus, in the long term we expect a situation where there is either a large amount of criminal 
activity or a complete lack of it. 

On the other hand, the system is monostable if the parameters are chosen so that 

/(0) = /(l) = 0, /'(O) >0, /(u) >0, for^ze (0,1). (13) 

For the monostable case we have that (0,0) is unstable, as lims_^o+ Q^5'(0) > 1) and {l,g{l)) is 
stable. The case when there is a single steady state solution is the least interesting and we do 
not consider it here. 

We first prove that there exists a traveling wave solution of ([s]) that connects two steady 
states: a hotspot invading zero-criminal activity region. Let us recall [8j that for the single 
parabolic equation ([T]) there exists a unique monotone traveling wave front solution S{x — ct) 
to ([T]), which satisfies 

S"{z) + cS'{z) + f{S) = 0, lim S{z) = 1, lim S{z) = 0. 
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Furthermore, S'{z) < for \z\ < oo and c > iff 

[' f{s) ds > 0. (14) 
Jo 

The monostable case differs as traveUng wave solutions exists for a semi-infinite interval of 
speeds. In particular, there exists a c* such that there are traveling wave solutions U{x — ct) for 
any c > c*. Here, it is always the case that steady state U = 1 invades C/ = 0. 

We have the following comparison principle. 

Lemma 1 (Comparison Principle). Let L > 0, (si,ui) and {s2,U2) be two solutions of ^ such 
that si(x,0) > S2(x,0) and ui{x,0) > U2{x,0). Then si{x,t) > S2{x,t), ui{x,t) > U2{x,t) for 
all t > 0. 

Proof. The function g{s) satisfies g{si) — g{s2) > /3(si — S2), with /3 as in ([2]), thus, s = si — S2 
and u = ui — U2 satisfy 

St > As — s + aL{x)u 
ut> fis — u. 

Since aL{x),(3 > the maximum principle implies that s{x,t) > and u{x,t) > for all 
t > 0. □ 

Finally, we recall the definition of supersolutions and subsolutions used here, for functions 
with discontinuous derivatives. 



Definition 1 (Supersolution). A function ^+{x) G C^(]R\ {0, L}) is a supersolution to (10) if it 
satisfies dxx''P+ + fi'>P+) < and for ^ £ {0, L} 

lim tp'^{x) < lim Tp'_^{x). (15) 

A subsolution, ■iIj^{x), is defined similarly by reversing the inequalities above. The above 
definition is easily extended to a finite number of discontinuities. 

2.2 Summary of Results 

In this section we summarize our results. Our main interest is the case when a population is 
said to have anti- criminal tendencies, which mathematically corresponds to the case the system 
([s]) is bistable and we discuss this case first. 

2.2.1 Bistable System 

We first study the case of a homogeneous medium (i.e L = 0). We show that there exist traveling 
wave solutions, {ip, (j), c) with c G M"^, such that for all 2 G M 

' ip"{z) + c^'iz) - ^{z) + a(t){z) = 
c</)'(z) + 9(V'(^))-0(^) = O . . 

O<V(^),0(^)<1 ^ ^ 

V'(-oo) = 1, (/.(-oo) = (?(1), ^(+00) = 0, (^(+00) = 0. 

Theorem 1 (Traveling Wave Solution for the Bistable System). Let L = and f{s) = —s+ag{s) 



up to translations, to il6h, Furthermore, 



be bistable, that is f{s) satisfies (12). Then there exists a unique {ilj{z),(f){z),c*) with c* G M, 



(i) c* > iff (14) is satisfied 
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(ii) ip'iz) < and 4>'{z) < for all \z\ < oo. 
Now, consider initial data that roughly has the shape of a traveling wave. In particular, 

liminf so(-2) > a, limsupso(^) < a (17a) 
liminf ^0(2;) > ^(a), limsupuo{z) < g{a), (17b) 



for z G M. The next result states that initial conditions that satisfy (17) lead to the eventual 
invasion of criminal activity everywhere. 

Theorem 2 (Homogeneous Case: Invasion of Criminal Activity and Stability). Let L = and 

{s(x,t),u{x,t)) be solutions to ([S]) with initial conditions, {so{x), uo{x)), such that < uo{x) < 
g{l) and < So{x) < 1 satisfying (17). Furthermore, let {ip{z), (p{z), c) be the traveling wave 
solutions. Then there exists zi,Z2 G M and bounded, positive functions qu{t) and Qsit) that 
approach zero as t ^ 00 such that 

i'{z + Z2) -qs{t) < s{z + ct,t) < ip{z + z*) + qs{t) (18a) 
(/){z + Z2) - Quit) < u{z + ct, t) < (l){z + z*) + Quit), (18b) 

for all t > 0. If additionally, the initial conditions satisfy: 

\uo{z) - (p{z + zq)\ < e and \so{z) - ^l^{z + zo)\ < e (19) 

for some zq, e > 0. Then there exists a function, u}{e), which satisfies lime_s.ow(e) = and 

\u{z + ct,t) — (j){z + zo)\ < u){e) and \s{z + ct,t) — ip{z + zq)\ < uj{e), (20) 

for all t > 0. In particular, 

lim u(x, t) = g{l). 



The inequality (20) tells us that once a solution is close to a traveling wave solution it will 
remain close for all time. Numerical simulations indicate that the solutions of the initial value 
problem converge exponentially fast to traveling waves. The next result is for the heterogeneous 
problem ((fsl) with L > 0). 



Theorem 3 (Unique Entire Solutions). There exists a unique entire solution (up to a time 
shift) {s{x,t),u{x,t)) to ^ defined onR xR, such that < s{x,t) < 1, < u{x,t) < g{l), 
st{x, t) > 0, and ut{x, t) > for all {x, i) G M x M. Furthermore, if {ip{z),<j){z), c) are the unique 
solutions from Theorem^ then 

s{x.,t) ^ il){x — ct)., u{x,t) ^ (p{x — ct), (21) 

as t ^ —00 uniformly in x £M. 

Such solutions are called transition fronts, and one interpretation for the prevention of the 
invasion of criminal activity is that the transition waves tend to a steady state and do not 
propagate as i — )• +00. 

The comparison principle for ([s]) enables us to reduce the problem of blocking the generalized 
traveling fronts, {s{x,t),u{x,t)), to finding a steady state solution, sl{x), that is monotonically 
decreasing. In particular, we require that the steady state solution satisfy 

lim sl{x) = 1 lim s'l{x) = 0. (22) 
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(a) (b) 

Figure 1: Two simulations with the same parameters but with different values of L, both where 
ran with a mesh-size of Ax = .01. On the right (6) shows the first numerical value that blocks 
the wave, L* = .85, on the left (a) the simulation was ran with L = .84. 



The condition at negative infinity allows for the initial propagation of the wave and the condition 
at positive infinity guarantees that the transition fronts are blocked. Thus, the problem is 

s'-fL{x,s) = Q (23) 

where, 

„ , . f f'(s) X < and x > L 



with the boundary condition (22). 

Our next result, and the crux of this work, is on the existence of critical gap length L*, such 
that for L > L*, transition fronts are blocked. On the other hand, for L < L* the fronts are 



delayed, but eventually propagate (see Figure la). Before stating this result, we define h G (0, 1) 
such that F{h) = 0, where 



F{s) := r f{v)dv. 
Jo 



Theorem 4 (Blocking the Generalized Traveling Fronts). Let {s{x,t),u{x,t)) be solution to ^ 
from Theorem\^ then there exist a unique L* > such that 

(i) if L > L* then the propagation of the waves are blocked, i.e there exists a steady state 



solution sl{x) to ([8|) that satisfies (22) 



(ii) if L < L* then the waves propagates. That is, for any e > there exists x^ > L and t^ > 0, 
such that 

s{x,t) > 1 — e for X > Xe, t > t^ (25a) 
u{x,t) > 1 — e for X > Xe, t> t^. (25b) 



Furthermore, L* is characterized as the unique value such that (23) with L = L* has a unique 



monotone decreasing solution, sl*{x), satisfying sl*{L*) = b and s'^*(L*) = 0. 

The results of Theorem [4] are in agreement with numerical simulations, see Figure [T| which 
shows the evolution of two numerical solutions with the only difference being the lengths of the 
gap. On the left, L = .84 and t = 350, one observes that after some time the wave propagates. 
Furthermore, the traveling wave front is eventually reformed after some time. On the right, 
L = .85 and t = 800 and, contrary to the previous case, the invasion of criminal activity is 
prevented. 
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Finally, we address the question of splitting resources for the bistable system. Given that 
the minimum amount of resources required to stop the propagation of the generalized wave is 
L = L*, we now ask the following question: if we split the resources into two regions with decay, 
that is two intervals with lengths Li and L2, which are separated by a distance d, does there 
exist a split and a distance d between these two intervals such that Li + L2 < L* is sufficient 
to stop the propagation of the waves? See Figure [2] for an illustration of the problem. The next 
result states that splitting resources can never help minimize the total amount of resources used. 

Proposition 1 (Splitting Resources). Let L < L* , with L* from Theorem Consider an 
arbitrary split of an interval of length L into two sub-intervals of length Li and L2 (as de- 
scribed above). Then, for any distance d > between the two intervals of length Li and L2 the 
generalized traveling waves will propagate. 



2.2.2 Monostable System 

In this section we discuss the case when the population's natural tendencies towards crime is 
indifference (this corresponds to so = 0). Mathematically, this corresponds to a monostable 
system. We first state the corresponding result for traveling wave solutions in the monostable 
case. 

Theorem 5 (Traveling Wave Solution for the Monostable System). Let L = and f{s) be 



monostable, there exists a c* £ M"*" and solutions {ipc{z) , (pdz)) with z = x — ct to (16) for any 
c> c* . Furthermore, when c < c* , such waves do not exist. 

An interesting conclusion is that assuming that sq < was crucial in preventing the prop- 
agation of crime. In fact, when the population is neutral towards criminal activity, then no 
amount of limited resources can prevent the propagation of crime waves. 

Theorem 6 (Monostable Criminal Activity). Let L > and (s(x,t), u{x,t)) be solutions to 
([s]) with f{s) monostable, then the waves propagate. Ln particular, (25) holds for some x^ > L 
and tf > 0. 



2.2.3 Nerve-Pulse Propagation Equation 

As noted earlier, system ([s]) is related to equation ([T]). In fact, our methods and results carry 
over to the bistable ([T]) , and the results of [T3] are included in Theorem |4j We can also qualify 
the invasion process by proving a similar result to that in Theorem |3j adding to the picture of 
the gap problem for the single parabolic bistable equation. 

Theorem 7 (Unique Entire Solution). Let L > 0, there exists a unique entire solution s{x,t) 
to ^ onRxR with < s{x, t) < I, st{x, t) > for all (x, t) G M x M. Furthermore, 

s{x,t) -i/jix + ct) ^0, (26) 

as t ^ —00. 



3 Bistable Case: Ant i- Criminal Activity Tendencies 

This section is devoted to proving the results for the bistable case. 



3.1 On Traveling Wave Solutions for the Homogeneous Problem 

Proof. (Theore m [I| The existence and monotonicity of the traveling wave solution, {(j),il),c), 
which satisfies (16), follows from the general result in [20] (Theorem 2.1) on the existence of 
traveling waves for monotonic systems. It is straightforward to verify that the conditions of this 
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theorem are satisfied. To prove part (i) we let S{z) be the travehng wave solution of ([7]) with 
unique speed c. Recall that it satisfies 

S" + cS' + f{S) = 0. (27) 

Consider the moving coordinates z = x — ct and define 'w{z,t) = s{z + ct,t) and v{z,t) = 
u{z + ct, t). Then any solution to v = {w, v) to satisfies A^v = with 

f Wt — Wzz — cWz + w — av 
Nv = < , ' (28) 

[ vt - cvz - g{w) + v. 

Let c > and assume for contradiction that c < 0. Applying the operator N to 
VI = {S{z),g{S{z))) gives 



iVv 



1 



-5"(z) - cS'iz) - f{S{z)) 
-cg'iS{z))S'{z)-giSiz))+giSiz)). 



Given that S'{z) < and that S{z) satisfies (27) we see that -S"-f{s) < 0. Also, -cS'{z) < 0. 



Thus A^vi < 0, which implies that vi is a subsolution. We just need to check that there are 
not issues at infinity. To do this, recall that the traveling waves are unique modulo translations. 
We take z* large enough so that 

S{z + z*) < ijj{z) and g{S{z + z*)) < (p(z) for z < R, (29) 

for some real number R and 



S{z + z*) < 6, for z> R, 



for some small 5. 



We now show that (29) actually holds for all z £R. RecaU that c > and S' (z) , ip' {z) , (p' (z) < 



and that we assumed that c < 0. Thus, 



Thus, w = S — ip satisfies 



with 



r{z)+fWz))<o 

S"{z) + f{S{z)) > 0. 



w"{z) + c{x)w > 0, 



5- V 

Note that c{x) < if w > 0, as S < 6. We apply the strong maximum principle to the set 
{w > 0} and conclude that w < since w = on the boundary, which is a contradiction. Thus, 
u; < and so S{z) is below 'iIj{z) for all z SM. A similar argument works for g{S{z)) and (j){z). 
Finally, since S and g{S) are traveling waves moving to the right and c < implies that the 
traveling wave solutions {iIj,4>) are moving to the left, this leads to a contradiction. The same 
argument works when c < and this proves (i). □ 

3.2 On the Dynamics of the Bistable Homogeneous System 

This section is devoted to proving Theorem [2j The proof of this theorem is a generalization 
of the method of [8J and the idea is to produce a suitable supersolution and subsolution in the 
moving coordinate frame z = x — ct. 
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Proof. (Theorem [2]) Let {^^{z), (l){z), c) be as in Theorem [T] and chose q^, q^, z* such that 
hminf so(a;) > 1 — > a hminfiio(2;) > g{l) —qu> g{o), 

and 

^{z + z*) -qs< sq{z) (j){z + z*) -qu< uq{z). 



This is possible by choosing z* positive and large enough as {sq{z) , uq{z)) satisfy (17). Define 
v_ = (u;_,t;_)"^ where 

W- = m.ax.{Q,il){z + £^{t)) — qs{t)} and u_ = max {0, + ^(t)) — . 

The functions ^{t),qs{t), and quit) will be defined later. Let t = z + ^(t) and obtain 



em'ir) - q'sit) - r{r) - ciP'ir) + V(r) - - a{<PiT) - qu{t)) 



iVv_ 



Making use of ( |16| we obtain 

/V^ _ / - - Qsit) + ag«(t) 

I e'(t)0'(T) -q'u + gWr)) - gWr) - - q.{t). 

Chose (5 > smaUso that for ■0 G [0, J]n[l-(5, 1] and (p £ [0,5]n[g{l)-6,g{l)] and < qs{t) < q^ 
then 

g{iP{T)) - gii^ir) - qs{t)) < Kqsit) 

with K satisfying 

OK < 1. 

Such K exists given the definition of g{s). If S,'{t) > in this regime 



(30) 



I'sit) - qsit) + aquit) 
Ut) + Kqsit)-quit). 



Solving the system 



Q'sit) = -qs{t) + aqu{t) 
I'uit) = -qu{t) + Kqs{t), 

with initial conditions ^^(O) = q^ and qs{0) = qg we obtain solutions 

qs{t) = Ae^-* + Be^+^ 

Quit) = 



(31a) 
(31b) 



a 



a 



where A-t = — 1 it ^/aK < and A, B satisfy 



1 

Jan 
a 



JOLK 

a 



1 





' A ■ 




' Qs ' 




B 







Inequality (30) guarantees that qsit) < and quit) < qu- From (31) we see that qsit), quit) > 0. 
Furthermore, 

lim qsit) = and lim quit) = 0. 
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We are left to treat the case when (-0, (p) G (5, 1 — 5) x ((5, g{l) — 5). Here ip', </>' < —7 for some 
7 > 0. Since j3 is the maximum value of g'{s), our goal is to find ^(t) such that 



iVv_ < 



-ieit)-q'sit)-qsit) + aqu{t) 
-li'{t)-q'^{t)+ pqs{t)-qu{t). 



By definition of qs{t) and qu{t) the first inequality above is automatically satisfied; therefore, we 
need only solve 

li\t) = -q'St)+Pqs{t)-qu{t) 
= {fi-K)qs{t). 

with initial condition ^(0) = z* . The solution is 

^{t) = zi + Aie^-* + Bie^+ 

where, Ai = ^^^-^'^'^ and Ai = ^^"^ zi = z* — {Ai +Bi). Note that £,{t) is increasing, but 

reaches a limit as t approaches infinity, 

lim ^{t) = Z2. 

Therefore, z* < (,{t) < Z2 and at t = we have 

■0(2; + z*) - < so(x) and 4>{z + z*) - qu < uo{x), 
which assures that v_ is a subsolution. Similarly, we can check that (/>(z+Z2)+'?u(i) > v{z,t) and 



'ijj{z + Z2) + qu{i) > w{z,t). Thus, we obtain the upper and lower bounds in (18). Furthermore, 



as the upper and lower bounds for u{x,t) in (18) approach g{l) as x,t ^ 00. Finally, to prove 



(20) we take ^^(t) \z* — Z2\ or oder 0{e). This is possible due to (19) and suffices to prove 



the result. □ 

3.3 Unique Entire Solutions for the Heterogeneous Problem 

This section is devoted to the proof of Theorem [3j which is a result for the heterogenous problem 
L > 0. Given v = (w, v)'^ a solution to ^ satisfies A^v = 0, where 

[ vt- g[w) + v. 

The proof of Theorem [s] requires some asymptotic estimates on which we state below. 

Lemma 2 (Asymptotic Estimates). Let {ip,(j),c) be defined as in Theorem^then the following 
estimats hold for and (j): 

76"^^ < 4^{z),(piz) < 6e-^' z>0 (33a) 
76^^ < 1 - V(^) < Se^"^ z<0 (33b) 
je^^' < gil) - giPiz)) < Sei^^ (33c) 



and for ip'^z) and (j)'{z) : 



-76"^^ < i^'{z),p'{z) < -6e-^' z>0 (34a) 
-76^^ < 'ip'{z),(t>'{z) < -lef"^ z<0, (34b) 



for 7, 7, 6, A positive constants. Furthermore, X> fi. 
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The proof of Lemma ([2| is standard so we only give the proof of A > //. (see [H [TT] for the 
derivation of the estimates for the single equation). We mention, however, that the conclusion 
that A > ^, which is not an automatic part of standard exponential decay estimates, will be 
very important later. 

Also, since the traveling wave solutions {ip{z), (piz)) are unique modulo translations, without 
loss of generality we set 

^(0) = I and 0(0) = ^. (35) 

The proof of Theorem [3] follow techniques of [1] . For the existence part of the proof it is useful 
to define, following [4J, the auxiliary function ^(t) by 

where c is the speed of the waves from Theorem [T| A is defined in Lemma [2j and M a positive 
constant to be chosen later. Note that ^(t) is well-defined on i G (— oo, — T) with T := ^^^^^'^^ ■ 
Furthermore, limt^_oo '^(i) = and 

C(t) =Me^("*+«). (37) 

For uniqueness it is useful to consider the spatial region where the front of the waves are located 
at a given time. To make this more precise, given G [Oi 5) we define the time function 

F^{t) := {x < : ry < s(x, t) <l-r], g{r]) < u{x, t) < g{l - r/)} . 

Since the waves are moving to the right there exists a time, Trj G M, such that Frf{t) C {x < — 1} 
for t G ( — 00, ?jj). We state an extension of Lemma 3.1 in [3] and leave out the proof as it follows 
similarly to that of Lemma 3.1. 

Lemma 3. For any r] G [0, 1/2), there exists a (5^ > such that for {s,u), the entire solutions 
to then 

st,ut > forx G Tr,{t), t G (-00, T^). (38) 

We now prove Theorem [3} 

Proof. We break the proof into four steps. First, we show that the piecewise function, v-|- = 
(u'+,u+), defined by 

w^(xt)-( 'f'^'-''^^'^^'-^ "^^^ andt;^fxt)-| + < ° (39) 

where z+ = x — ct — ^(t) and z- = —x — ct — £,{t), with a suitable ^(t) (which approaches zero 
as t — >• —00), is a supersolution for t G (— 00,— T), for T large enough. Note that w^{x,t) is 
independent of the spatial variable in the gap and is composed of a pair of self-annihilating 
fronts. Second, we find a suitable time range where v_ = {w-, v-) defined by 

w.{x,t) = lfy+^-'^^y-^ .ndv.{x,t) = l t^y^^-'^^y-^ (4o) 

^^[0 x>0, ^^[0 x>0, ^ ^ 

for y+ = X — ct + ^(t) and y- = —x — ct + ^{t), is a subsolution. In the third step, we take the 
v_ at a suitable negative time as initial data and construct a monotone sequence of solutions, 
which are bounded below by the the subsolution and bounded above by the supersolution for 
values of t that are negative enough. From this we conclude the existence of an entire solution, 
that is bounded above by u+ and below by V- for negative enough values of t. By the definition 
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of the super and subsolutions then we obtain (21). The final step to show uniqueness and the 
proof will be combination of the uniqueness proof found in [3] for the entire solution and the 
proof for Theorem [2j 

Step 1 (Supersolution): We aim to show that is a supersolution in a suitable time range. 
Case 1: Consider a; < and apply M to v-|_ 



m = ! -eitw{z+)+^'{z^)) 

\ + (p'iz^)) + gWz+)) + g{^{z-)) - g{^j{z+) + V'(z-)). 

Since "tp^z) is always negative the first term in the above equality is strictly positive. For the 
second term we make use of the inequality |llj . 

\g{a)+g{h)-g{a + b)\<Lah, (41) 

for any < a, 6 < g{l) and some constant L > 0, obtaining 



If 2;+ < then |x| > —ct — ^(t) as we are considering x < 0; thus, z_ = —x — ci — ^ > 0. Thus, 
in this case 2+ < and z+ > 0. We now invoke Lemma [2] we obtain 

Since X > fj, the above quantity is positive provided 

M7 > L7. (42) 
On the other hand, when 2+ > (note also that this implies that z„ > 0) 

-C'{t){^'{z+) + 0'(z-)) - L^{z+)ijiz-) > M7e^('=*+«)e-^(^-^*-«) - L-f^e-^^''-''*-^^e-^'^-^-'^-^^ 

Hence, for x < if 

M7 > L72, (43) 

we have that Mv-|_ > 0. 
Case 2: If x > then 



-2(c + e{t))^'{-ct - e(t)) + 2'iPi-ct - C{t)) - 2a^{-ct - ^(t)) 

-2(c + e{t))<p'{-ct - m) - g{m-ct - m)) + 20(-ct - m) 
-2em'{-ct-m) 

-2i'{t)4>'{-ct - m) + 2g{^{-ct - m)) - g{2ij{-ct - m))- 



Let Ti > be sufficiently large so that t G (-00, -Ti) then -ct-^{t) > 0, then 7p{-ct-^) < f. 
Since, g{2il){—ct — ^(t))) = in this case Thus, v+ is a supersolution for t G (— 00,— Ti) as 
g{s) < for all s < sq. 
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Step 2 (Suhsolution): The case when x > is trivial thus assume that x < 0. Applying 
M to v_ 



i-c + eW{y+) - - 9Wy+) - V'(y-)) + Hy+) - Hv-) 

^ / eWiy+) - ^P'iy-)) 

\ i'{(p'{y+) - (t>'{y-)) + g{^{y+)) - g{^{y-)) - g{^{y+) - ^{y-))- 

Observe that y+ < as x < and y_ > because x < 0. There are two possibilities we need 
to consider, y+ < and y+ > 0. In the latter case, both ip{y+) and ipiy-) are bounded above 
by ^. First, this implies that ~ i^'iy-) ^ and 4>'{y+) — (t>'{y-) < 0. Second, it implies 

that — g{xp{y-)) — g{ip{y-f-) — ■ip{y-)) = 0; hence, in this case Mv_ < 0. We are left to 

consider the case when y+ < 0, which we work out next (once again using (41)) 

^ ~ [ Me^'^^*+^\-Se^y+ +je-^y-) + LiP{y_){iP{y+) - ip{y-)) := //. 
Rewrite the first term 

/ = _Me^(^+''*+^) ^^gM(-ct+C)+x(/i-A) _ ~g-A(-ct+c)' 

Therefore, if t < — Ti, replace Ti by a larger value if necessary, then 

l^^^{-ct+m) _ ^^-K~ct+m) > for t G (-oo, -Ti] (44) 

because A > /i. For the second term we use the bound 

LV(y-)(V(y+)-V'(y-))<^<5e-^^-, 

which gives, 

// < -Me^(^-+^*+«)(5e^(-^*+«) -7e-^(-^*+«) -L(^e-^(-^*+«)) 
<0, 

provided Ti is chosen negative enough to satisfy 

^eM(-rf+e) _ ^e-A(-rf+0 _ i5e-A(-ci+0 > q fort e (-oo, -Ti] 

(with Ti replaced by its updated value). 

Step 3 (Entire Solution): Let n be a positive integer and {sn{x,t),Un{x,t)) be the solutions 
to ([s]) for X G (—00,00) and t G (— n,oo) with initial data 

Sn{x, —n) = 'W-{x, —n), Un{x, —n) = f_(x, — n). 



From the definitions of V-|- , v_ , see ( 39 ) and ( 40 ) , respectively, we see that 

w^{x, —n) < Sn{x, —n) < w+{x, —n) and f_(x, —n) < Vn{x, —n) < w+(x, —n) 
By the comparison principle 

u;_(x, t) < Sn{x, t) < w+(x, t), v-{x, t) < Un{x, t) < v+{x, t) for (x, t) G M X (— n, — Ti). (45) 
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In particular, (45) hold for t = —n + 1 = — (n — 1), which implies that for x G M 

s„_i(x, -n + 1) = w-{x, -n + 1) < Sn{x, -n + 1), 
u„_i(x, -n + 1) = -n + 1) < u„(x, -n + 1). 

Another application of the comparison principle gives that 

Sn-i{x, t) < Sn{x, t) and Un~i{x, t) < Un{x, t) for (x, t) S M X (— n + 1, — Ti). 

Thus, we obtain monotone increasing sequences and which enables us take 

the limit as n — t- oo and obtain solutions {s{x,t),u{x,t)) to ([S]) on {x,t) G M x M. Furthermore, 
we have the additional bounds 

W-{x, t) < s{x, t) < W-^-{x, t), v-{x, t) < u{x, t) < u+(x, t) for {x, G M x (— oo, — Ti). 



The definitions of v_|_, v+ also give (21). We are left to prove that 

st{x,t),utix,t)>0. (46) 

We begin by looking at the time derivative of v_ 

dtw. = {-c + ^'){i;'{y+)-i^'{y.)). 

Take t negative enough so that — c + ^'{t) < and arguing as before {ip'{y+) — ip'{y-)) < 0. 
Note that if y+ > this inequality is clear and if y+ < then we take t negative enough). Thus, 
dtw^ > for t sufficiently negative. Using (45) we obtain that dtSn{x,t), dtUn{x,t) > if n is 
large enough and thus the maximum principle gives that 

dtSn{x, t) > 0, dtUn{x, t) > for t G {—n, oo). 



Taking the limit as n — >• oo and using the fact that dts{x, — oo) > (same for dtu) proves (46). 
Step 4 (Uniqueness): Assume for contradiction {s{x,t),u{x,t)) is another solution to (Isl) that 



also satisfies ( |21[ ). Since (s, u) and (s, u) both satisfy (21 ) then for any e > there exists a time 
G M such that 

||s(-,t) -s(-,i)lloo < e ||u(-,t) -n(-,i)||oo < e (47) 

for all t G (— oo,te). For to £ {—co,t^ — e] define 

5+(x,t) ■.= s{x,t + to + eC{t))+eqs{t), U+{x,t) := u{x,t + to + eC{t)) + equ{t) 
S~{x,t) := s{x,t + to - €^{t)) - eqs{t), U~{x,t) := u{x,t + to - e^(t)) - equ{t). 

As in the proof of Theorem [2] we aim find ^(t), qu{t) and qs{t) such that {S~^,U~^) are 
supersolutions and {S~,U~) are subsolutions for t G [0, — to — e]. Since the proof is very 
similar to that of Theorem [2] we only summarize the steps for the subsolution and note that the 
supersolution is checked similarly. Let t = t + to + €^{t) and apply the operator Ai (defined in 

(Q) to {s-,u-y. 

dtS- - + S-- aU- = -e{e{t)dts{T,x) + q',{t) + q,{t) - aqu{t)) 
dtU- - g{S-) + U~ = -ei'{t)dtu{T,x) - eq'^{t) + g{s{T)) - g{s{T) - eq,{t)) - eq^it). 

In the region outside the front of the wave, that is ifx ^ J-",y(t + to + eC(i)) then for tj small enough 
we have that g{s{T)) — g{s{T) — eqs{t)) < Kqs{t), where once again we chose k such that 

QK < 1. 
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Then, as in the proof of Theorem [2] we can find positive, bounded functions qu{t) and qs{t) 
that satisfy (31). Recall, that both qs{t),qu{t) — )• as t — )• oo. On the other hand, when 



X € Tr]{t + to — e^{t)) then we invoke ( |38| ) noting that this holds when t <T^ — — e, and solve 
for i{t) 



at) 



13 -K 



sit) m = o- 



Here, since /3 > k we also get ^(t) is increases but approaches a limit as t — )• oo, say 

lim ^(t) = M. 

Thus, we get that indeed {S~,U~) is subsolution and {S^,U^) is a supersolution for t G 



[0, — to — Me]. Inequality (47) implies that 

S'{x,0) < si{x,to) < S^{x,0) and {x,0) < ui{x,to) < U^{x,0). 

Implying that for t £ [0, — to — e] then 

S~{x,t) < si{x,t + to) < and U'{x,t) < ui{x,t + to) < U^{x,t) 

for all X G M. By a change of variables t = t + to, and dropping the tilde, this is equivalent to 

s(x, t - e^(t - to)) - msit - to) < s{x, t) < s{x, t + e^(t - to)) + eqs{t - to) 
u{x,t- eC{t - to)) - equ{t - to) < u{x,t) < u{x,t + e^(t - to)) + eg'„(t - to), 

for t G [to,T^ — Me] and to G (— oo,T^ — Me]. Taking the limit to — )■ — oo, using the facts that 
qu{t - to),qs{t - to) and ^(t) M, then 

s(x, t - Me) < s{x, t) < s{x, t + Me) 
u{x, t- Me) < u{x, t) < u{x, t + Me) 

for (x,t) G M X (— oo,rr, — Me]. Finally, letting e — )• gives that s = s and u = u and we 
conclude. 

□ 

3.4 Preventing the Propagation of Criminal Activity 

This section is devoted to the problem of preventing the invasion of the criminal activity. Con- 



sider the alternative representation of the steady state solution. Any solution to (23) satisfies, 
for any yi,y2 G (-oo,0] or yi,y2 G [L,oo), 

1 ry2 ry2 rsL{y2) 

7, [(s'Liyi))' - {s'L{y2))'] = - I 4ix)s'dx) dx = I f{sL{x))s'dx) dx = / /(0) d9. 



yi 



yi 



Letting yi = — oo and 1/2 = above then (22) implies 



1 



is'dx)f 



f{9) dO. 



for any x < 0. On the other hand, when x G [0, L] the solution must have the form 

s^(x) = yle~^ + Be" 
for some ^4, i? G M. As we seek solutions the following must hold: 
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\{B-Af = T{A + B)-T{l) 
sl{x) = Ae-^ + 5e^ 



— oo < X < 
< X < L 



(48) 



F{Ae~^ + Be^) - F{sl{oo)) L < x < oo. 



If there are values L G R"^ and A,B gM. with the bound < A + i? < 1, then one can exphcitly 



build a solution to (23) using the quadrature method (see [7J) for appropriate values of sl(oo). 



Indeed, if iI^l{x) is the unique solution to 

u" + f{u)=Q xg(-oo,0) 
ii(_oo) = 1 u{<d) = A + B. 

Such a solution exists as A + B = s(0) € [0, 1] and it is monotonically decreasing; in fact, 



iI)'l{x) 



12 1 f{e)de. 



So tp'i^ix) < because the right hand side of the above equation is bounded away from zero. 
Thus, if a unique solution, 7pR{x), exists to 



u" + f{u) =0 X £ (L, cxd) 

u{L) = Ae~^ + Be^ n(oo) = sl(oo). 



(49a) 
(49b) 



then we can build an explicit solution of the form 



iPl{x) X G (— oo, 0] 

sl{x) = { ^e-^ + Se^ x£ {0,L) 
iPr{x) X G {L,oo]. 



Note that (49) has solutions for certain values of sl(oo) 



Remark 1. To find solutions to (23 ) it suffices to find appropriate values of L, A, B and a unique 



iPr{x) satisfying (49). A blocking solution will necessarily have that sl(oo) is bounded away 
from one. 

Now we move to the proof of Theorem [4| In order to prevent the propagation of the waves 
the existence of a steady state solution that is uniformly bounded away from one as x — )• oo is 
necessary. In fact, we will show that with sufficient resources a strictly decreasing solution that 
approaches zero as x — t- oo will always block the propagation of the wave. We first show that 
for L sufficiently large such a steady state solution exist. 



3.4.1 Case I: Sufficient Resources 



Our goal in this section is to build super and subsolutions to (23) for L large enough, which 
both decay to zero at positive infinity. In particular, we seek a supersolution, (/)+(x), such that 
4>+{Q) = 1 and (p+{L) = b, where b is defined by F{b) = 0. To be precise, let 4>+{x) = Ae~^+Be^ 
for X G [0, L], we seek to find A, B, and L such that 



A + B = l, Ae~^ + Be^ = 6, and - Ae'^ + Be^ 



0. 



This reduces to finding L such that coshL = |, which always has a solution because < 6 < 1. 
The unique values are given by 

Lo = cosh"^ f M >0, B= ^e-^\ and A = ^e^° . 



We use Lq to construct a monotonically decreasing solution to ( 23 ) with L = Lq 
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Lemma 4 (Base Case). There exists a monotonically decreasing solution, slo{x), to ^ with 
L = Lq, such that 

lim SLn{x) = 1 lim slo(x) = 0. 

X— >— oo a;— ^-oo 

Proof. We prove this lemma in three steps. The first step is to construct a suitable supersolu- 
tion and subsolution. We then build a solutions by using the supersolution and subsolution as 
barriers. The first solution constructed is not necessarily monotone, but it is bounded above by 
b for all X > L. Using this upper bound we prove there also exists monotone solution that is 
below the first solution we constructed. 



Step 1: (Super /subsolution) Let ^/'(x) be the unique solution to Cauchy problem 

-u" = f{u) X G [Lo,oo) 
n(Lo) = 6, u'{Lq) = 0. 



Such a solution exists and in fact following the arguments from the beginning of ^3.4 we can 
show that '4)'{x) < and lim^j-^oo i^ix) = 0. Now, define 4'+{x) by 

{1 X < 

b (giog-^ + e-^oe"") < X < Lo 
Tp{x) X > Lq. 

The function (/>+(x) is a supersolution as it satisfies 

lim (j)'^{x) > lim <f>'_^_{x), 

a-— a;— >5+ 

for ^ = and ^ = Lq. Similarly, we build a subsolution, 0-(x). Let tp- now be the solution to 
the Cauchy problem 

-u" = f{u) xG[-oo,0) 
u{-oo) = 1, n(0) = 0. 

As before, we know that such a solution exists, and we define the subsolution 

ip- (x) X < 
X > 0. 



Furthermore, s'(0) = —y"^ Jq f{^)d 0, which is bounded away from zero and so the condition 
on the derivative at x = for a subsolution is satisfied. 



Step 2: (Barrier Method) Having constructed a supersolution and a subsolution we use the 
standard barrier method (see for example [18j) to show that there exists a solution to (23). This 
steps follow those of the proof existence of the entire solutions in Theorem [s] (see Step 3 of the 
proof), so we leave out the details. 



Step 3: If the solution is monotone then we are done. If it is not monotone then the lack of 
monotonicity happens in the gap. Indeed, for x > L we know that sl(x) < b and then 

-lis'dx))' = l^'^^\{9) de<o. 

Since, the solution must be concave up in the gap then it reaches a minimum at a unique 
point X* G (0, L). If sl(x*) > a we construct a supersolution by letting s+(x) = sl{x) for x < x*. 
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s+(x) = sl{x*) for X £ [x*,xi], where xi is the unique xi > L such that sl{xi) = sl{x*). For 
X > xi then s+{x) = sl{x). Hence, Sj^{x) cuts off the non- monotone part of the solution sl{x). 
The other possibihty is that sl{x*) < a and in this case cut off from with x < L and 

y > L being the unique points with sl{x) = siiu) = a. In both cases, s+{x) is a supersolution. 
Once more the barrier method will allow us to construct a solution slo{x), which is monotone. 
Indeed, in the latter case, sl{x) < a for x > x* and s^^(x) < in this regime. In the former 
case, SLg{x) cannot have a non- negative derivative everywhere because it can only be zero (and 
thus later negative) if slq{x) = b, which can never happen. This concludes the proof. □ 

The maximum value of the solution at at L is indeed 6, which is the maximum value that 
the homoclinic orbit through the point (s(x) = 0, s'{x) = 0) attains. For L = Lq is is clear that 
the steady state solution has a value less than or equal to b at Lq as the solution must remain 
below the supersolution. This brings up the question of wether there exists an L < Lq such that 
s{L) = b and s'{L) = 0, so that the maximum value can actually be achieved? The answer to 
this question is provided in the following lemma. 

Lemma 5 (Necessary Resources). There exists a unique L* > such that there exists a mono- 



tone solution, sl*{x), to (23) with sl*{L*) = b and s'j^i,{L*) = 0. Furthermore, any solution, 



sl{x), to (23), with L < L* , must satisfy sl{L) > sl*{L*) = b. 



Proof. The first objective is to find L such that sl{L) = b and s'^{L) = so that we know that 
Ae~^ + Be^ = b and Ae~^ = Be^. Note that this immediately satisfies the last equality in (48) 
with sl(oo) = 0. Substituting this into the first equality gives 

hbsinhLf = J^{b cosh L) - J-(l), 

which has a solution. Indeed, for L = the right hand side is positive and the left hand side is 
zero. On the other hand, for Lq = cosh~^(j) the left hand side is positive and the right hand 
side is zero. Hence, by continuity there exists a solution, which we call L* G (0,Lo). Given L* 
then we have a solution 

s\{x) = ^(e^*-^' + e^-^*) for x G [O,^'^], 

which is monotonically decreasing and reaches a minimum at x = L*. Also, for any x > L* we 
have 



\ {d^s\x)f = / f{e) de, 



b 

(x) 



which implies that s*{x) < b for all x > L*. Furthermore, since dxS*{x) ^ when s G (0,6), 
s*(x) cannot be a constant on any interval, and dxS*{x) is continuous then dxS*{x) < for 
all a:: > L with s*(x) ^ 0. On the left side of the interval, s|^(0) > b and so for x G (— oo,0] 
sl(x) G (6, 1) and s'{x) < 0. We are left to prove that 

lim s^(x) = 1 and lim s'^{x) = 0. 

a-— >— oo a;— ^+oo 

This is true, of course, only because s'^{x) > b for x < and < 6 for x > L*. 

Claim: lim2;_!._oo s^ix) = 1. We know that the solution is bounded by one and the it 
is monotonic. Thus, we only prove that limit is not somewhere in between b and 1. This is 
accomplished by proving that the reaction term f{u) approaches zero. Consider the solution 



w{x) 



19 



which is the solution to 

= 1 in (-1,1) v = Oon{-l,l}. 

For X < 0, define 

r{x)=mf{f{e):ee[b,sl{x)]}. 

We win prove that r(x) < ^jjit^Top"' '^^i^h then implies that /(s) vanishes as x — )• — oo. Assume 
for contradiction that this is not the case, so there exists a xq G (— oo, 0) such that 

Note that xq cannot be a local minimum, as dxxS'^ixo) = — f {sl{xo)) < 0, thus we can find xi 
close to Xq such that s^(xi) < s^(xo) and r < dist(xi, 0) < dist(xo, 0). Here, r is chosen so that 
(50) still holds, that is. 

Now consider the solution. 



r(xo) > ^. 



z{x) = r{xo)r'^v 



X — Xl 



to 



— z" = T{xo) for X G Ir{xi) 
z = on 91, 

where Ir{xi) is the interval centered on xi with radius r. We now scale z{x) by r by defining 

= Tz{x). 

It is clear that for r small enough z'^{x) < sl*{x) for x G It{xi) and as r increases there is a 
first r* and x* where z'^ touches the graph of s^(x), that is sl*{x*) = z'^ (x*). We know also 
that X* is inside the interval. Thus, using the fact that the maximum of z{x) is achieved at xi, 
we obtain 

SL^ix'') = Z-^\x*) < T*z{xi) = l!£^^ < SL.{xi) < Sl*{xo) < 1- 

Thus, r* < '^s 2 < 1, by our assumption on r. Define w{0) := t*z{x) — s*r (x) on /^(xi). From 
the above argument we know that 

w <0 e Ir{xi) and w{x*) = 0. 

Now, since sl*(x*) < sl*{xq) there exists a neighborhood of x*, V, such that s|^(x) < s|^(xo) 
for all X G V. This implies 

s'L*{x) < -r(xo) for X G V. 

So in this neighborhood, V, 

w"{x) > r*(xo)(l -r) > 0. 

However, this contradicts the fact that w{x) has a local maximum at x*. Thus, r(x) — )• so 
f{s) —7- oo and from this we conclude the claim. The remaining limit, lim^_^_|_oo s^(x) = is 
proved similarly by looking at the equation that 1 — s|^(x) solves. 

Now, consider the case when L < L* clearly s*(x) is a subsolution to the problem 

s" + /l(s,x) = 0. (51) 



Any solution to (51), sl{x), must satisfy sl{x) > s*{x); hence, it is true that sl{L) > s*{L) > 
s*{L*) and in particular sl{L) > b. Similarly, for L > L*, s*(x) is a supersolution to (23) and 
thus sl{L) < s*{L) < s*{L*). Hence, L* is unique. □ 
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Next, we show that for any L larger than L* there also exists a monotone decreasing solution. 

Lemma 6 (Sufficient Resources). There exists a monotonically decreasing solution, sl{x), to 
(23) when L > L* , such that 

lim sl(x) = 1 lim sl(x) = 0. 

Proof. Having constructed a solution for L* , call it s*{x) we now show that for L > L* there 
also exists a steady state solution. Indeed, note that we have 

Slx + fL*{x,S*)=Q, 

and so s* is a supersolution to the problem, 

s" + fLix,s) = 

as fiix, s*) < fL*{x, s*). The rest of the proof follows the proof of Lemma |4j □ 
3.4.2 Case II: Insufficient Resources 

Now, we treat to the case when L is small and the wave propagates. Indeed, when L = we 
know that for appropriate initial conditions the wave will propagate, see Theorem [TJ In the first 
lemma we prove the existence of a symmetric solution, connecting one to itself, for any L > 0. 
In the case when L = the symmetric solution is s{x) = 1. 

Lemma 7 (Symmetric Steady State Solution). For all L > there exists a solution, Smix), to 
(10b), with the following properties: 

(a) Sm{x) is symmetric about x = L/2. 

(h) lim^^ioo Sm{x) = 1. 

(c) Sm{x) is monotone increasing on [L,oo). 

Proof. (Lemma [T]) Since this problem is translation invariant, we prove the existence of a sym- 
metric solution by translating the problem to the left by L/2. In other words. 



a{x) 



a X £ (-00, -2^] n [f ,oo) 



x£ {-^ ^ 



2 

2 ' 2 )' 



SO that any symmetric solution must be of the form Sm{x) = Ae~^ + Ae^ on [-^, As 
noted before the steady state problem is equivalent to (48) with J^{v{oo)) = J-"(l), and so the 
symmetric condition implies that the following equality must hold 



2 \ K'iJ J V V 2 

Letting z = 2^cosh(^) the above reduces to 



^^tanh^ =F{z)-F{l), 



which always has a solution as —J-{\) > 0. Hence, we can build explicitly the symmetric solution 
about the origin by the method discussed earlier (see Remark [l]) such that 



lim Sm{x) 

x— >±oo 
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We are left to show that this solution is monotone outside of [— -L/2,L/2]. To do this we prove 
that s'jy^{x) > for x > L/2 given the fact that lim-r^oo Sm{x) = 1 (by symmetry then < 

for X < —L/2). Chose 6 such that Sm{L/2) < 1 — 6, then the limit implies that for (5 > there 
is a constant R > 0, sufficiently large, such that 

Smix) > I — S for X > R. 

Consider a translated solution 

s'^{x) = Sm{x + r) for r > 0. 
Note that s'^{x) > 1 — 5 for x > R — t, so for t = R — L/2 we obtain 

sl,{L/2) > s^{L/2). 

Invoking the maximum principle for semi-infinite domains for bistable reaction-diffusion systems, 
see [2], we obtain that 

sj^(x) > Sm{x) for X > L/2. 

Now, define 

T* = inf {r > I slj^ix) > Sm{x) for x > L/2} , 
it is clear that r* < R — L/2. Let us assume for contradiction that r* > 0, by definition of r* 

G := inf(s;;(x) - > 0, 4 = [L/2,R]. 

If G > then, by standard elliptic theory [10], there exists a e > such that 

sl^{x) > Sm{x) xela, te {t* - €, T*). 

The maximum principle for unbounded domains on [R,oo) gives that for some eo G (0,e), 

s;^''°{x) > s^{x) X e [L/2,oo). 
This contradicts the definition of r*. If G = 0, there exists a sequence {^jtlfcgpj S la such that 



We translate the solutions by defining sj^(x) = Sm{x + Ck), and obtain a bounded sequence of 
solutions. Thus, a subsequence converges to a limit, Smix), on compact sets. Furthermore, as 
G = we have that (0) = Sm(0) and sjj^(x) > Sm{x) for all k; thus, Sm{x + T*) = Sm{x) for all 
X > by the strong maximum principle. Thus, Smix) must be a periodic function with period 
r*. However, recall that G [L/2, R] and so s^{x) — )• 1 as x — )• oo, which gives a contradiction. 
Thus, r* = and we conclude that Sm{x + 1) > Sm{x) for all t > 0. 

□ 

We now show that for small L there are no solutions that are bounded away from one for 
X > L. 



Lemma 8. For L < L* there are no non-negative solutions, sl{x), to (10b) that are bounded 
above by one, and such that 

lim sl{x) = 1 lim sl{x) < 1 — e. (52) 

X— oo x~^oo 

for any e > small. 
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Proof. (Lemma Is]) Assume, for contradiction, that there exists a solution sl{x) to (10b) such 



that satisfies (52). Since L < L* from Lemma [s] we know that si{L) > b. Furthermore, for 
L < y < X 

{My))' - {M^))' = 2 r^^^ f{0) de. (53) 

Note that the integral is bounded away from zero as sl{x) G (6, 1); thus, s']^{x) / s'j^iy)- Hence, 
the derivative is monotone. Assume for contradiction that s'j^{x) < for all x > L, then 



sl{x) < sl{L), from (53) we get that s'j^{L) < s'^(x). Thus, limx^oo s'j^ix) = 0, if sl{x) is to 



stay above b. Looking x — )• oo, we obtain 



My))' = 2 lim f{e)de<o. 



This is a contradiction, and so it must be that s'{x) > and by prior arguments, see the proof 
of Theorem [5} we know that 

lim = 1. 

This gives the final contradiction and we conclude. □ 
3.5 Proof of Theorem S] 

Now, we have all the tools necessary to prove our main results. We begin with the proof of 
Theorem |4l 

Proof. (Proof of Theorem |4| Let s{x, t) and u{x, t) be the unique solutions to ([3| with initial 
data that satisfies (17). Since s{x) = 1 and u{x) = g{l) are supersolutions for ^ and (0,0) are 



subsolutions, then by standard parabolic theory s{x,t) — )• s{x) and u{x,t) — t- u{x) uniformly on 



compact sets, as well as its derivatives, where s{x) and u{x) solve (10). Furthermore, given the 
initial data st{x,t), ut{x,t) > 0. Therefore, s{x) / and u{x) ^ 0. First assume that L > L* 
from Lemma |4] and Lemma [5] we know that there exist monotonically decreasing steady state 
solutions 

lim sl{x) = lim g{sL{x)) = 0. 

This proves (i). Now, consider the case when L < L* . In this case, the only steady state solution 
is the symmetric solution. In fact, from Lemma |8] we know that all steady state solutions must 
satisfy sl{x) £ (^,1) with s'{x) > for x > L. This automatically implies that the steady 
state solution for the moving average of crime, call it ul{x) also satisfies ul{x) £ {g{b),g{l)) 
with u'j^{x) > for X > L. [Start with initial conditions such that st{x,t),ut{x,t) > 0]=^ the 



inequalities (25) are satisfied. In fact, we have that 



lim s{x,t) = Sm{x) lim u{x,t) = g{smix)). 



□ 



From proof above we immediately obtain as a corollary the same result for equation ([T]) with 
bistable /(s). Recall, that this result was conjectured in [13] . 

Corollary 1. Let L* be as in Theorem^ and s{x,t) a solution to ([7]) with f{s) bistable then 
for the following hold 

(i) If L > L* then the wave propagation is blocked. 

(ii) If L < L* then the wave propagates. 
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3.6 Is Splitting Resources Useful? 



In this section we focus on the issue of weather sphtting resources can improve the situation by 
minimizing the total resources needed to prevent the propagation. To make this rigorous let us 
set L < L* and arbitrarily split L into two regions of lengths Li and L2 separated by an interval 
of length d. To summarize we have 

Li + L2 := L < L*. 

Now, the reaction term, f{u), is replaced by decay term, — s, in the two regions of length Li 
and L2. The question we aim to answer is this: Is there a distance d between the two gaps of 
length Li and L2 so that the wave propagation is prevented? For clarity, we state the problem 
mathematically by defining a new function. 

f(^,^-jfis) (-oo,0)n(Li,Li + d)n(L + d,oo) 

Jd[x,s)-<^_^ xG [0,Li]n[Li+d,L + d], ^ 



and solve the new problem 



s" + fd{s) = 0. (55) 



Figure [2] gives an illustration of the double gap problem. It is self-evident that if the distance 
between the gaps is very large then, by definition of the critical length, the wave will propagate 
through the first gap. If the traveling front passes and is allowed to reform, once it hits the 
second gap the length will once again be too small to stop the propagation of the wave. Hence, 
we know that for d too large splitting the resources make will necessarily increase the total 
amount of resources required to prevent the propagation of the waves. 

Figure 2: Double Gap Problem 
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Proposition 2. Let {s{x),u(x)) be a solutions to (55) with Li + L2 = L < L* , then for any 
e > there exists an x^ > L + d such that s(x) > 1 — e and u{x) > g{l) — e for all x > x^. 

This proposition implies that splitting resources will never improve the situation, or reduce 
the total number of resources. In the first step we prove the result for small d. The second step 
is an inductive step, where prove if the result holds that for an arbitrary d > 0, then it holds 
for d + e for e > 0. 

Proof. We prove this in two steps. 

Step 1: Assume that d < L*—L and let s*{x) be the strictly decreasing solution obtain in Lemma 
[5} In particular, s*{x) satisfies s*^ + fL*{s*,x) = 0. In this case we have fL*{s*,x) < fd{s*,x) 
which implies that s* is a subsolution of 

s" + fd{x,s) = 0. 

Note in particular that the inequality is strict as s*{x) £ (b,!) for x G [0,i^*]. Therefore, 
s{x) > s*{x) and u{x) > g{s*{x)). In particular, s{L*) > b which implies, as we saw before, that 
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the wave propagated and (25) is satisfied. 



Step 1: (Inductive Step) Next, assume that given di > arbitrary the traveUng waves propa- 
gate. It suffices to show that separating the gaps by a distance of di + ei, for ei > 0, will not 
help. We will prove this with the help of suitable subsolutions that we now define. Let Sd^{x) 
be the solution to (55) with d = di, which by assumption satisfies (25) for some > L + di. 
Furthermore, let Ss{x) be the symmetric solution to 

-s" = -s L + di < X < L + di + ei 
s{L + di) = Sd^L + di) s'{L + di) = s'^^{L + di), 



and iIj{x) is the solution to 



-s" = f{s) x> L + di + ei 
s{L + di + ei) = Sd{L + di). 



Both Ss{x) and ip{x) exist and furthermore il){L + di + ei) > b by definition of S(i^(x). Then 

s{x) = 



Sdiix) x<L + di 
Ss{x) L + di < X < L + di + ei 
iIj{x) X > L + di + ei, 



is a subsolution to (55) with d = di + ei. Indeed, we have that for x < L + di 

-s" - fdi+ei{s,x) < 0, 
and for X > L + di we have that fdi+ti {s, x) = fd^+ei {s, x) 

s" - fdi+eAs,x) = 0. 

Hence, s{x) is a subsolution and given that s{L + di + ei) > b then any solution, s{x) must 
also satisfy s{L + di + ei) > b. Since, the splitting was arbitrary then previous arguments this 
concludes the proof. 

□ 



4 The Monostable Case: Neutral Criminal Activity 

Theorem [5] directly from Theorem 2.2 in [20]. In terms of preventing the invasion of criminal 
activity, we prove that for all L > all steady state solutions increase to one as x, t — )• 00. In 
fact, in this case, due the strength of the instability of zero the solutions will always approach 
the symmetric solution unless the initial data is exactly zero. 

Proof. (Proof of Theorem |6]) Consider the initial value problem 

s" + /(s) = xG[L,oo) 
q{L) = (3, 

for any (3 £ (0,1]. All bounded solutions approach one as x ^ 00. As lims{x,t) = s{x) and 
limti(x,t) = g{s{x)) on compact sets then it must be the case that linix^t^oo s{x,t) = 1 and 
limt^x^oo u{x,t) = 5(1). 

□ 

As before we obtain the exact result for monostable single parabolic equation. 

Corollary 2. Let L > be arbitrary and s{x,t) a solution to ^ with f{s) monostable then the 
wave always propagates. 
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5 The Single Parabolic Equation 



The proof of Theorem [T] follows from Theorem 2.1 in fl] if we study the symmetric problem, 

St - Sxx- fLix,s) = (56) 

where, 

. , . f f(s) X < -L and x > 
fLi-,s) = ['_\^ -L<x<0. 

Additionally, one needs to assume that the invasion is from the right. Hence, the traveling wave 
solution we are interested in satisfies for z = x + ct with c > 

r r{z)-ci^'{z) + f{z) = o 

\ V(-oo) =0 V(+oo) = 1. ^ ' 

Note that this problem is equivalent (symmetric) to the previous problem. With the help of ip{z) 
in the pre-invasion process we show that there is an entire solution on M x M, which converge 
to ip{z) as t — 7- — cxD. As the proof is similar to that of Theorem ^ and follows the proof of 
Theorem 2.1 in [4J we leave the details to the reader. 



6 Discussion 



We have studies the effect that a population's natural tendencies towards crime have on criminal 
activity patterns. The simplest case is when there is a natural tendency towards criminal activity, 
equivalently in our model Sb > 0. In this case, in the long run one expects either a constant 
hotspot or a warm-spot. In the case when there is indifference, Sf, = 0, lack of criminal activity 
is a steady state; however, it is unstable and if the payoff is high enough, criminal activity will 
have a tendency to dominate. The most interesting case is when the population has a natural 
tendency to be peaceful, or to avoid criminal activity. Here, there is an interesting interplay 
between the natural tendency and the payoff of committing a crime. In this situation, there are 
three steady states, complete lack of crime, a small amount of criminal activity, or a hotspot. 
The lack of criminal activity and hotspots steady state solutions are stable and the small amount 
of crime is unstable. Hence, in this case it will almost always be either all or nothing, so that in 
the long run there is a high criminal activity or zero criminal activity. 

We focused on the neutral and anti-crime tendencies to study the invasion of criminal ac- 
tivity into areas that start out with zero criminal activity. In both cases, we prove that the 
invasion of criminal activities is possible, under the right circumstances. This invasion happens 
via traveling wave solutions. Note that this also holds for M" in general, of particular interest 
in this application is the case. Indeed in Figure [s] we see a traveling wave solution for s{x,t) 
moving in the x-direction. In the neutral case, or monostable case as we have been referring 
to it, the high criminal activity steady state solution is always more stable and the invasion 
will always be that direction. On the other hand, in the anti-criminal activity case, or bistable 
case, it is easier to have the high crime steady state solution be more stable; however, there is 
the possibility that the invasion to be in the other direction. More specifically this implies that 
there are parameters that lead to zero crime areas invading the high crime rate zones. 

We have also addressed the issue of the prevention of the propagation of criminal activity 
by applying a finite amount of resources to reduce the payoff of committing a crime to zero on a 
finite interval. This corresponds to employing crime deterring strategies such a providing more 
police officers to the area, decreasing opportunities, etc. In the bistable case, we showed that 
there is a critical length of the interval with zero payoff coefficient, or equivalently a minimum 
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Figure 3: 2D Traveling Wave of the Propensity to Commit a Crime 
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amount of resources, required to prevent the propagation of criminal activity. The amount of 
resources required depends on 7^(1), vaguely measures the stability of the hotspot steady state, 
as was shown numerically. An interesting conclusion of our work is that in the monostable case 
there is no possible amount of finite resources that will stop the propagation of crime. Hence, 
to stop the propagation of crime it is essential that the population be naturally anti-criminal 
activity or to start out with zero-crime. 

There are many interesting questions left to answer. For example, in the anti-criminal activity 
scenario there are two ways to prevent the invasion of the criminal activity waves. The first is 
to have coefficients such that / f{s) ds < and the second is to employ enough resources on an 
area of length L* to diminish the payoff to zero there. An interesting question is to determine 
which of these options is more cost effective. In the former case corresponds to changing both 
the payoff to commit a crime and the natural crime-tendencies, this could be done say via 
preventive educational efforts. The latter is simply a direct intervention in reducing the payoff 
of committing a crime. 



A Lemma [2] 

We have shown that there exists traveling wave solutions, {ip{z), 4>{z)), with speed c > which 
satisfy 



ip"{z) + cil)'{z) +TP + a<j){z) = 
c(t)'{z)+g{xl,{z))-(t>{z)=Q 



(59a) 
(59b) 



This can be turned into a system of ODE's by defining (v,q,w) = (V', </>, V'O- Then the above 
system is equivalent to 



V = w 



9{v) 



w = —cw + V — aq. 
The Jacobian at a steady state v* = {v* ,q*,w*) is 



J(v*) 



1 

-a'jv*) 1 Q 

c c 

1 —a c 



with eigenvalues 



X-'-) (A^ + cA-1)-^^. 
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Recall that ^'(0) = and so for the steady state (0,0) the eigenvalues are 

, 1 , -c ± Vc^ + 4 
Ai=^, A± = 2 ■ 

It is useful to define the functional 

Q(A, k) := (^\ - (A^ + cA - 1) - A;. (60) 

As k = —Si^ > then we see that indeed A > /x. 
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